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Abstract 



a. 

I , In this review we present a biased review of the ground state properties of the 

Falicov-Kimball models in 1, 2 and oo dimensions, considering either fermions or hard 
core bosons. In particular we want to show the very rich structure that these models 
£j ' exhibit, and to point out the analogies and differences associated with the statistic of 

the quantum particles and the nature of the lattice (bipartite or not). The flux phase 
problem is also discussed. 
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1 Introduction 



1.1 Definition of the model 

The Falicov-Kimball model is a lattice model of quantum particules interacting with clas- 
sical particles. 

The lattice A = {x} is defined by a finite or infinite set of sites x in M d . Usually one 
considers d-dimensional cubic lattices Z d , and A C However it is interesting to consider 
more general lattices, e.g. the 2-d triangular lattice which, contrary to cubic lattices, is 
not bipartite. Indeed, particles on bipartite lattices have particle-hole symmetries, and 
one would like to know what properties remain valid for systems which do not have these 
symmetries. As we shall review, the rigourous results which have been obtained so far 
concern the cases where the coupling between quantum and classical particles is either 
very strong (in any dimension), or very weak (for d=l). For intermediate coupling, the 
only results concern reduced phase diagrams (sec. 1.4) and have been obtained using 
exact numerical methods. To investigate this intermediate coupling regime, it appears of 
interest to study first the model on the Bethe lattice in the limit where the coordination 
number becomes infinite; indeed in this limit one obtains a mean field model which should 
be simpler to study rigorously. 

In the original model, the quantum particles, called "electrons" in the following, are 
spinless fermions described by creation and annihilation operators a x ,a x , satisfying the 
usual anticommutation relations 

{a x ,a y } = {<4, a*} = 0, {<4, a y } = 5 x>y . (1) 

It is well known that the Pauli principle is important to get a crystalline state, since such 
a state can not occur if fermions are replaced by bosons Q. However hard-core bosons 
do indeed lead to crystaline structures Q and they are interesting to study. Spinless 
hard-core bosons are described by operators a x ,a x satisfying the following relations 

{al,a x } = 1, {4)4} = i a x, a x} = ( 2 ) 

but 

[a*,a*] = if x^y {a* = a) x ,a x ) (3) 

For the finite system A, the kinetic energy of the quantum particles is defined by the 
operator 

a;,j/£A 

where t xy = t* x = \t xy \e t6xv . Complex hopping constants t xy are introduced to model 
particles in an external magnetic field: given a circuit C = (x\,X2,-- - ,x n ,x n+ i = x\), 
i.e. an ordered sequence of sites Xi in A such that t Xi)Xi+1 7^ for all i = 1, . . . , n, the 
magnetic flux through this circuit is defined by 

n 

= 9 =<>i,Xi+l mod 27r - ( 5 ) 
i=l 
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In the original Falicov-Kimball model, only hopping between nearest-neighbour sites were 
considered, with 



(6E iflx-^l 
10 otherwise. 

Since then, several extensions of the model have been investigated, introducing quantum 
particles with "spin", described by operators a xa ,a xcr , together with the kinetic energy 

k a = - e E< a -V'- ( 7 ) 

x,y£A <j,<j' 

For example, in the standard Hubbard model 

K A = ~Y^ ^txyaLaya- ( 8 ) 

x,y£A <r=T,J. 

For the asymmetric Hubbard model 

Ka = ~ E E ^«Lv ( 9 ) 

x,yeA <r=1,l 

with the Falicov-Kimball model given by 

t X y = txyi ^ X y = 0- (10) 

For the Montorsi-Rasetti model 

^a = -E E ^2/4 CT v ( n ) 

x,yeAa,<r'=1,l 

and it has been shown that this model is equivalent to the Falicov-Kimball model. 

The classical particles, called "ions" in the following, are assumed to have a hard-core. 
They are described by random variables W x G {0, l},x £ A, where W x = 1 means that 
the site x is occupied by an ion and W x = means that the site x is not occupied by an 
ion. In the spin-language, the random variable is s x , where s x = 2W X — 1 € { — 1,+1}. 
The classical particles could be ions, impurities, spins, localised f-electrons, depending 
on the physical system under investigation. 

Again several extension have been considered. For example in the static Holstein 
model s x € IR (and is associated with phonons), while in the static Kondo model s x = 
(s x ,s x ,s x ) G IR, |s x | = 1, describes localized magnetic impurities. 

The classical particles do not move and they do not have any kinetic energy. However 
in the definition of equilibrium states one will take annealed averages over all possible 
configurations of ions. 

Except for the hard-core between ions, and the statistics for electrons, it is assumed 
that the only interactions are between electrons and ions when they occupy the same 
lattice site. It models for example the screened Coulomb interaction U which can be either 
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positive (repulsion), ore negative (attraction). In other words, given the configurations 
W = {W X } X £A, or s = {s x } X £\, of the ions, the interaction is described by the operator 

H^(W) =2UJ2 (W x - ±) (a\a x - ^ (12) 



or, equivalently, 



H^(s) = uY,s x Ua x -^\. (13) 



The factor 1/2 is introduced for convenience to exhibit symmetry properties. In the canon- 
ical formalism it amounts simply to adding a constant; in the grand canonical formalism 
to a redefinition of the chemical potentials. In conclusion, given a configuration of the 
ions, the hamiltonian of the finite system A is given by 

H A (s) = K A + H^(s) (14) 

to which is added in the grand canonical formalism the contribution 

- Ve^ala x - (j,i-^2(s x + l), (15) 

xeA xGA 

with /i e ,/ij the chemical potentials for electrons and ions. 

Again several extensions can be considered. For example, for electrons with spins the 
following hamiltonian has been considered 



H A (s) = K A + uY J Sx(n x -^\+U'Y J KT a *T " \) ( a U a *l ~ \) 
xeA ^ ' xeA ^ J Vi " J 



(16) 



where n x = o^a xT + a^a^. 



1.2 Interest of the Falicov-Kimball model 

One of the most fundamental problem in condensed matter physics is to understand the 
phenomenon of phase transitions, in particular why all elements and many compounds 
crystallize in periodic structures. It is well known that one of the driving principles behind 
this ordering is associated with the structure of the ground states; however it is still not 
clear what physical mechanisms are necessary for phase transitions to occur, especially 
for quantum systems. Indeed, because of the quantum fluctuations, the ground state 
properties, as well as the low temperature behavior, are difficult to extract for quantum 
systems. 

In 1969 Falicov and Kimball || introduced their model to study metal-insulator transi- 
tions in mixed valence compounds of rare earth and transition metal oxides as an effect of 
the interaction between localised f-electrons ( "classical particles" ) and itinerant d-electrons 
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(quantum particles). Later this same model was considered to study ordering in mixed 
valence systems, order-disorder transitions in binary alloys, and itinerant magnetism. 

This model was reinvented in 1986 by Kennedy and Lieb ||] as a primitive model for 
matter to study crystallisation. In this interpretation ions are represented by classical 
particles and itinerant electrons by quantum particles. Depending on whether there is 
one, or several, electronic bands near the Fermi level, one is lead to consider neutral 
(equal density of ions and electrons), or non-neutral systems. For p e = pi = 0, which is 
the symmetry point of the system and corresponds to the neutral case with p e = pi = i, 
Kennedy and Lieb proved the existence of long range order (crystal) at low teperatures for 
any coupling U and any dimension d > 2, together with the absence of such order (fluid) 
at high temperatures. Their result was then extenden by J. L. Lebowitz and N. Macris in 
1994 ||, for values of chemical potential in a neighborhood of the symmetry point. 



1.3 Problems 

The first step to study phase transitions is to obtain the zero temperature phase diagram, 
i.e. the ground states. 

Problem 1 : "Ground states" (T = 0) 

Find the configuration of ions which minimize the energy. 

Two approaches have been used to solve this problem. 

In the canonical formalism, given (N e ,Ni), the number of electrons and ions, the 
problem is to find the configurations of ions s = {s x } which minimize the energy, i.e. 

E Ne (s)= min A E Ne (s') (17) 

s>: Ni(s')=Ni 

where iVj(s) = i ^2 x ^\(s x + 1) and Ej^ c (s') is the infimum of the spectrum of the hamil- 
tonian H^^ e (s') restricted to the N e electrons sector. 

In the grand canonical formalism, given (p e ,pi), the electron and ion chemical poten- 
tials, one starts from the partition function. Introducing an effective free energy T by 
means of a partial trace over the electrons degrees of freedom, i.e. 

Z A (p,p e ,pi) =J2Tr{e-^ HA ^- fMeNe -^ iNi ^} = J^ e -^A(«;^e,w) ( 18 ) 

s s 

one is lead to study a classical spin lattice system. Defining the effective hamiltonian for 
this classical spin system by 

Ea(s; fi e , Pi) = lim J 7 A(s;l3,pe, Pi) (19) 

the problem is to find those configurations of ions s which minimize E\(s'; p e ,pi), i.e. to 
find the zero temperature phase diagram of the effective hamiltonian. 

Since the interaction is on-site, the quantum-mechanical problem for fermions can be 
solved by looking at the eigenvalues ej(s) of the hamiltonian h\(s) for one electron moving 
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in the potential U (s x ) defined by the configuration of ions. Therefore for fermions systems 
E A (si; n e , m) = ^ i e j( s ) ~ A*e] - y ^2 s x ~ ^iNi(s) (20) 

ej(3)<fi e x 

= E Ne (s)-n e N e - f i i N i (s) with iV e = N e (n e ). (21) 

On the other hand for systems of bosons with hard cores, we have a truly many-body 
problem of interacting particles; in this case one has to proceed via the effective free energy 
and the limit (3 — > oo. 

Problem 2 : "Flux phases" (T = 0) 

Let $ = {0c}) where <f>c denotes the magnetic flux through the elementary circuit C. 
The problems are the following: 

i) Given {p e , pi,Q), find the configurations of ions s which minimize H\(s;&, p, e , /jj). 

ii) Given (// e , /jj), find the configurations s and the fluxes $ which minimize 
H A (s; $,{i e ,fj,i). 

Problem 3 : "Low themperature phase diagram" 

The real problem in the study of phase transitions is to show that the zero temperature 
phase diagram is stable at low temperatures. 

This question is discussed in the talks by R. Kotecky, D. Ueltschi, and N. Datta, and will 
not be considered in this review (see ref. [7]- [13]). 

1.4 Methods 

Several methods have been introduced to study the zero temperature phase diagram. 

a) Expansion of E\(s; [i e , /Xj) in powers of \U\ 

This method has been introduced in both the canonical and the grand canonical for- 
malism. It is valid in any dimension, but is restricted to 

\U\ > ct and \i e e] - \U\ + ct, \U\ - ct[ (22) 

where c is some constant depending on the lattice and t = max \ t xy \. The conditions ( [22] ) 
imply in particular that this technique is restricted to 

"neutral" systems, i.e. p e = pi, if U < (23) 
"half-filled" systems, i.e. p e + p t = l, if U > (24) 

where = m and p% = are the (average) electrons and ions densities. 

In the case of fermions, the expansion of -Ejv e (s) m powers of \lf\~ 1 is easily obtained 
using (20), together with the following property: with z the maximal coordination number 
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of the lattice and \ U\ > zt, then, for any ions configuration s, the spectrum of the 1-electron 
hamiltonian h\(s) has a gap containing the interval ] — \U\ + zt, \U\ — zt[; for p e inside 
this gap the number of eigenvalues ej(s) < p e is iVj(s) if U < and |A| — Ni(s) if U > 0. 

Therefore, for iV e = iVi(s) if U < 0, and for N e = |A| - ATj(s) if U > 0, i.e. for /i e in 
the above gap, we have 

£N , W = 4g % + 5 L£ d ^{_|- M } (25) 

where C is a contour in the complex plane enclosing all negative eigenvalues. 

Iterating the resolvent identiy for [z — /ia(s)] -1 we obtain explicitely the desired ex- 
pansion 

E Ne (s) = T-^-E n (s) (26) 



n>l 1 1 



with 



n+l 



rp f \ _ V- (~l) m ("-1)1 TT I r97 N 

nlj " 2^ m ( m -i)i(n-m)! 1 ^ 1 ' 

x 1: ... ,x n+1 eA v 7 v y i=i 

where the sequence (x±, . . . , x n +i) must contain at least one emply site and one occupied 
site, and m = m(s) is the number of sites x%,i = 1, . . . , n+ 1, such that s x . = —1 if U > 0, 
and = +1 if 17 < 0. 

For hardcore bosons, the \lf\~ 1 expansion is obtained using the closed loop expansion 
of Messager-Miracle [12] for the effective free energy T{J5\ s, p, e , /Xj). Taking the limit 
j3 — > oo yields the expansion for E\(s, p e , p,j) which is convergent if the conditions fl22]) 
are satisfied, and implies the restriction Eq. (23), and (24). 

In conclusion for both statistics we are able to write the effective hamiltonian as 

E A {s; p e , m) = H^(s;fi e ,fii) + R^(s;n e ,m) (28) 

where 

k 

H^(s; Me , (a) = H<® (a; p e , + £ jTTj^W- ( 29 ) 

n>l ' ' 

At this point, the strategy is to study the ground states of the truncated hamiltonian 
and to control the rest R^ k \ With this technique one can prove that the phase 
diagram of the effective hamiltonian E\(s; p, e , /ij) is rigorously given by the phase dia- 
gram of the truncated hamiltonian, except for domains of width (which can be 
explicitely estimated) centered on the boundaries of the phase diagram of H^ k \ At this 
order nothing can be said concerning (p, e , jXi) in these domains. Going to order (k + 1), 
the ground states of will give the ground states of Ea(s; p e , /Xj) except for domains 
of width \U\~( k+2 \ and so on. 

This method has been applied up to order k = 3 to construct the phase diagrams for 
d = 1 and 2, cubic and triangular lattices, both statistics, for systems with or without 
magnetic fields. However it is restricted by the conditions p e = pi if U < and p e + pi = 1 
if U > 0. 
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b) Expansion of E\(s) in powers of \U\ 

This method has been used only in the canonical formalism. Contrary to the first 
method, it is valid for any rational densities (p e , P i) of electrons and ions, but so far it is 
restricted to 1 dimension only. 

Assuming t xy = t / only if \x— y\ = 1, then, using Rayleigh-Schrodinger perturbation 
theory, one obtains for the ground state energy density (for U < 0) 

2t 1 \W I 2 

e{s;p e )=2p e -- S mTrp e -Up ePt + —^^U 2 ln\U\+0{U 2 ) (30) 
7T 47Tt sin TTp e 

where p e = p/q (p prime with respect to q) and 
1 9-1 

W q (s; Pe ) = -J2e- t2 ^ x s x . (31) 

Given p e = p/q, the problem is then to find the configurations s which minimize e(s; p e ). 

c) Reduced phase diagrams 

This method is valid for arbitrary coupling constant U and arbitrary densities (p e ,pi), 
resp. (p e ,pi), but gives only qualitative results. The idea is to select some restricted 
class of ions configurations, e.g. all periodic configurations with period < 16 together with 
mixtures of two such periodic configurations, and to search, by means of exact numerical 
computation, for the configurations in this class which minimize the energy. This yields 
the so called reduced phase diagram. It has been applied in 1 and 2 dimensions. 

From these numerical calculations, one observes that this approach yields results con- 
sistent with those rigorously established in the limit of large or small U. Moreover the 
reduced phase diagram appears rather stable, i.e. increasing the class of configurations con- 
sidered, the boundaries of the previous reduced phase diagram become domains of smaller 
and smaller width, where new ground states configurations appear, while the main part of 
the diagram is not modified. The situation is similar to the case discussed above, passing 
from the truncated hamiltonian to H^ k+1 \ 

d) Finite systems 

In 1 and 2 dimensions explicit numerical computations of the energies have been de- 
velopped to find the exact ground state configurations of some finite systems. In fact, it is 
this early approach which led to the conjecture of "molecule formation" and to the idea, 
later proved, that for small \U\ and small densities (e.g. p e = pi < \) the ground state is 
not periodic. 
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2 Results for 1-dimensional systems 



In 1 dimension, using the closed loop expansion for the effective free energy, one first 
concludes that hard core bosons are identical to fermions 0. 

In the following we consider only the attractive case U < 0. From particle-hole sym- 
metry one then obtains similar conclusion for U > 0. 

Using the canonical formalism, the following results were obtained. 

Theorem 1 (Strong coupling |]T4|l ) 

a. Let p e = pi = p/q, with p prime with respect to q, then for \U\ > U cr (q) the ground 
state is the most homogeneous periodic configuration with period q. The position 
of the ions on the sites (0, 1, . . . , q — 1) is given by W x = 1 for x = kj where kj is 
solution of the equation 

pkj = j mod q j = 0, 1, . . . ,p - 1. (32) 



b. Let pi = bp e ,b ^ 1, then for \U\ > U seg (b) the ground state is the segregated 
configuration where all ions clump together. 



Theorem 2 (Weak coupling [15]) 



Given p e = p/q with p prime with respect to q, and \U\ <C q, then 

a. for pi £] — the ground state is a mixture of two periodic configurations s' 

and s" , with p' i = p" = 2-±i (p' e = p" e = p/q), where the position of ions is given 



by eq. (||), with j = 0, 1, . . . ,p' - 1 (resp. j = 0, 1, . . . ,p'). 
b. for pi = Pi/q, Pi not necessarily prime with respect to q, 

1. if Pi € ]0.37, 0.63[ ; the ground state is the periodic configuration of period q, 
given by eq. (|32]) . 

2. if pi < 0.37, or pi > 0.63, the ground state is a mixture consisting of a periodic 
configuration, with period q and ion density p\ = p'i/q, solution of eq. (|32|) , 
together with empty (p" = 0) or full (p'( = 1) configuration (and p' e = p" e = p/q), 
except for a countable set of densities (p e ,pi), with pi £ [|,0.37] U [0.63, f], 
where the ground state is the periodic ground state, with period q, given by eq. 
( |3~2"| ) . This countable set is given by the solutions of an equation, and we have 
for example 

117 6 14 

P l -P e - 4' 3' 20' 17' 39 
1 6 10 14 

Pi - 2p e 

This property shows that the periodicity of the pure phase is fixed by the electron 
density: it is the smallest period necessary to open a gap at the Fermi level. Let us also 
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remark that the critical ion densities given above (i.e. p\ = 0.37 and 0.63) are approximate 
values. The exact values are p c and (1 — p c ), with p c solution of 2i:p c = t&nirp c . 

The results of theorem 2 show that there is a close analogy with the Peierls instability; 
however this analogy is valid only for pi 6 [p c , 1 — p c ] and can be seen as follows: For U = 
(p e = p/q, Pi = Pi/q fixed) any ions configuration is a ground state and the probability to 
find an ion at a given site is uniform, equals to pj. This is the "undistorted state", which 
is metallic (no gap). For U 7^ 0, sufficiently small, a particular configuration is selected 
which has period q. It corresponds to the "distorted state", which has a gap at the Fermi 
level and is insulating. For pi £ [p c , 1 — p c ] the ground state is a mixture of a metallic and 
insulating phase (except for exceptional densities) and such a situation does not occur in 
the standard theory of Peierls and Frohlich. 

To extend these investigations to arbitrary values of the coupling constant U, and to 
avoid difficulties associated with mixtures, it is more convenient to work in the grand 



canonical formalism. Reduced phase diagrams were obtained in [16] for several values of 
U, by means of exact numerical calculations. Some examples are illustrated on figure 1. 





\ 1). 


;\ D c b 






\W\0(U " 3 ) 








0(U J )\^ 


D. 






(a) Strong coupling (b) Weak coupling 

Figure 1 



Figure la) shows part of the chemical potential plane for strong coupling, i.e. U < It. 
D-, D + , D c b represent domains where the ground state is the empty (p. L = 0) , the full 
(pi = 1), and the chessboard configuration. The parallel stripes corresponds to domains 
with p e = Pi = const. 

Figure lb) shows part of the chemical potential plane for weak coupling, i.e. U € 
[—24,0]. The "vertical" stripes with same denominator q correspond to periodic ground 
states with fixed electron density p e = p/q and different ion densities p, t = Pi/q. 

Combining all the results obtained so far, either rigorously or by reduced phase dia- 
grams analysis, one is led to the following conclusion 

Conjecture Given the electron density p e = p/q, p prime with respect to q, and the ion 
density pi = Pi/q, Pi not necessarily prime with respect to q, if the ground state is periodic 
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then it has period q and the position of the ions is given by the solutions of the equation 
pkj = j mod q, j = 0, 1, . . . pi - 1. (33) 
In all other cases, the ground state is either 

1. a mixture of two periodic configurations with pi = Pi/q and pi = {pi + l)/q given by 
the solutions of eq. ( p3| ) 

2. a mixture of one periodic configuration solution of eq. (^), together with either the 
empty (pi = 0) or the full (pi = 1) configuration. 

3. the segregated configuration. 

Let us remark that eq. (^) is reminiscent of the circle map theorem, and appears in 
many situations. 

3 Results for two dimensional systems 

The ground state properties of the Falicov-Kimball model on the square lattice, in the 
limit of strong coupling, are discussed in the talk by T. Kennedy. 

Reduced phase diagram analysis for the square lattice and arbitrary values of the 
coupling U has been conducted by Watson and Lemanski [fT^ ], |l8|| . Their results show 
that the properties discussed in sec. 2 are not specific to one dimensional systems, but 
also appears (together with new properties) in two dimensions. 

In this section we want to exhibit the difference between fermions and hard-core bosons, 
between square and triangular lattices, and to discuss the flux phase problem. 

This analysis has been conducted, within the grand canonical formalism, in the limit 
of strong coupling, using the \IJ\~ 1 expansion discussed in sec. 1.4. It is thus restricted to 
neutral (if U < 0) or half-filled (if U > 0) systems (ref. §, ||). 

The phase diagrams for fermions and hard-core bosons in the presence of an homo- 
geneous magnetic field, defined by its flux <3? through elementary cells (in particular it 
could be zero), are represented on figures 2a and 2b for the square lattice, and on figures 
4a and 4b for the triangular lattice. These phase diagrams are those of the truncated 
hamiltonians; the exact phase diagrams are identical except for small domains around the 
boundary curves which have been explicitely evaluated M. 



In fig. 4 (a) and (b), r_ and r + are the empty (p = 0) and full (p = 1) configurations; 
T5 and f 5 are the periodic configurations with densities ^ and | . 

The analysis of the triangular lattice has been extended to order 3. At this order the 
phase diagramms exhibit periodic configurations with densities p = (0, 1,5,3,3,3,3,5,5,5) 
as well as (1 — p). We note that there exist two different structures with densities 4, and 
with density \. The interested reader should consult the original article 0. 
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(a) Phase diagram to order 3 for fermions (b) Phase diagram to order 3 for hard-core 

on the square lattice bosons on the square lattice 



Figure 2 
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Figure 3: Configurations Si, $2, S3 which appears in fig. 2 (a) and (b). Si, §2, S3 are obtained 
by particle- hole transformation. S-, S+, S c b are respectively the empty (p = 0), the full (p = 1), 
and the chessboard {p = \) configurations. 

It remains to discuss the flux phase problem. In particular one would like to know in 
what cases, if any, the magnetic flux will decrease the energy of the system, and what is 
the optimal magnetic flux to obtain the state of minimum energy, i.e. the ground state. 

For hard-core bosons it is easy to see that the optimal magnetic flux is always zero 
(diamagnetic inequality) 0. 

For fermions the situation is more subtle; it depends on the lattice and the densities. 
On the square lattice, one finds that for density \ f ne optimal flux is 4> = tt through each 
plaquette, and the ion configuration is the chessboard structure; for densities | and | the 
optimal flux is no longer uniform but periodic with period 3, and = or 7r (the ion 
configurations are the structures S3, S3 of figure 5); similarly for densities | and |, at 
order 3 the optimal flux is non uniform, but periodic with <j) = or arbitrary (and the 
configurations are S\,S\); for densities and 1, the fluxes are arbitrary. We note that 
at this order the configurations with densities \ and | do not appear (but they might 
appear at the next order however). We also should remark that on the square lattice the 
optimal flux is maximum {i.e. 4> = ir everywhere) for density p = \, which is the maximum 
possible density, because of the particle-hole symmetry. 

The situation is very different for the triangular lattice. One finds that the optimal 
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(a) Phase diagram to order 2 for fermions 
on the triangular lattice. 



(b) Phase diagram to order 2 for hard-core 
bosons on the triangular lattice. 



Figure 4 



flux is maximum (cf) = ir everywhere) for densities p = I, | ; however for density 
P = § ' 2> §' ^ ne optimal flux is not uniform, but periodic, with (f> = or tt. (Let us 
recall that on the triangular lattice the particle-hole transformation is not a symmetry). 
Similarly for densities p = y and |, the optimal flux is not uniform, but periodic, with 
4> = or arbitrary (this arbitrariness may be lifted at the next order); finally for densities 
and 1, the fluxes are arbitrary. 

The results for the flux phase problem are illustrated on figure 5. 

We conclude this discussion with the following theorem Q. 



Theorem 3 For fermions systems 

1. for the configurations s+ = {s x = +1} and s_ = {s x = —1}, the effective hamilto- 
nian is independant of the magnetic fluxes. 

2. for any configuration s ^ s+ and s_, there exists Uq(s) such that for U > Uo(s), 
the optimal fluxes (i.e. those wich minimize the energy) are 

(a) for the square lattice 

= h ^ 2 Sit6S ° occupied (34) 
I otherwise 

(b) for the triangular lattice 

in J tt if or 1 site of A are occupied 

Pa = \ , . ( 35 ) 
(J otherwise 
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Square Lattice 



S— S\ S3 
1 1 h 



S4 
— I— 



5*3 Si S + 

H 1 1 * Mi ~ Me 



p = 



1 1 

5 3 



farbl 







arbl ;\ 7r 



7T 



0\ / 



tt I \ farb 



arb 



Triangular Lattice 



r- n r 3 
1 1 



H 1 h 



^ = ™ [arbl 1 ^ 



TT 



7T I \ 7T 



Figure 5: Optimal magnetic fluxes, [arb] means arbitrary at order 3, pi = p e = p. 



4 Preliminary results for oo-dimension Bethe lattice system 

In this section, the system is defined on the Bethe lattice with coordination number z 
(fig. 6) and by the hamiltonian 



(36) 



with W x £ {0, 1} and 



7= z if l*-yl = 1 

otherwise 



(37) 



We shall then consider the limit z — > 00. The Falicov-Kimball model on this infinite 
dimensional Bethe lattice has attracted a considerable interest since 1989 and the reader 



should consult the references [22], [23] for more informations. 
The density of electrons on the site x is expressed as 



(38) 
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Figure 6: Bethe lattice. The shells I = 0,1,2, 3 are represented. 



with u> n = j|(2n + 1) the Matsubara frequencies. 
The strategy proceeds as follows. Using 

1. Grassmann variables 

2. the fact that H is quadratic in a x , a x 

3. the fact that there are no closed loop on Bethe lattice 
and introducing the complex variable 

in = lU n - H e 

we obtain the explicit expression 

(W x ) - 1 (W x ) 



G xx {u3 n ) 



in + I x (U n ) in + U + I X (Uj n ) 

where 

I x {u n ) = \ t 2 Gy} X M 

y:\x-y\=l 

(A \ x is the Bethe lattice with the site x deleted) and 

-i 



l + e 



Y x ((3,He) = p J2(Htn + 4k)] " HZn + U + I x {uJ n )]} 
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4. At this point one takes the thermodynamic limit in a symmetrical manner, by in- 
creasing to infinity the number of shells. In this limit 

G«W = G|K), (W x )= ai (44) 

where I = 0, 1, 2, . . . denotes the shells (I = is the center point, I = 1 are the z 
nearest neighbours, and so on). 

5. Taking the limit of infinite dimension, z — > oo, we have from ( fll|) and (|40| ) 

/ x (w n ) = //(cj n ) = t 2 G i+ iK) (45) 

Gl{uJn) = e„ + * 2 WiK) " 6, + ^ + ^+xK) (46) 



«,= [l + e-fH/H-r-W,!*))] 1 
6. We can then take the zero temperature limit (3 — > oo, to obtain 



(47) 



1 

Yi(He) = 7r dLo{ln[£ + t 2 G l+1 (u;)}-lnlZ + U + t 2 G l+1 (uj)}} (48) 

7. We thus conclude from (47) that in the limit — > oo 

either a/ = and this happens iff Yi(fi e ) > fa (49) 
or ai = 1 and this happens iff Yi(fj, e ) < fa (50) 

8. The problem is reduced to the following ones: 
Given a = [ckq, ai,... , a p -i], a p G {0, 1}, 

(a) show that eq. (46) has a unique solution, Gi(u n ),l = 0, 1, . . . , n — 1 



(b) from Gi{uj n ) and eq. (45) determines Yi(fj, e ) 

(c) Using eq. (49) find those (fj, e ,fa) for which the "periodic" configuration a 
(periodic with respect to successive shells) is the ground state. 



With this procedure, the following result was obtained ([|20|, |2lJ). 
Theorem 4 For fi e e]2t, U - 2t[, then 

• for fa — fi e < —t 2 /U, the ground state is the empty configuration (W x = 0). 

• for fa — fa > t 2 /U, the ground state is the full configuration (W x = 1). 

• for fa — fa; € [—A, A], with 

t 2 t 4 t 6 t 8 

A =u-*m +K u' + ^ (51) 

the ground states are the period 2 configurations ex = [0, 1], and a = [1, 0]. 
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• for ^ — fi e G [-B, —A], and for fii — fi e £ [A, B] with 

t 2 i 6 t 8 

the ground states are the period 3 configurations given respectively by a = [0,0,1], 
[0,1,0], [1,0,0] and by a =[0,1,1],[1, 0,1], [1,1,0]. 

5 Conclusion 

Although extremly simple in its formulation, the Falicov-Kimball models exhibit a very 
rich structure in all dimension. In particular, on can prove 

• Periodic ground states 

• Molecules formation 

• Segregation 

• Coexistence of phases 

• Metal-insulator transition 

• Peierls instability 

• Devil staircase (maybe complete) 

• Farey tree properties 

• Flux phases 

Among the open problems, on can mention the following 

• Prove the conjecture for 1-dimensional systems 

• Find the exact equation for the segregated configuration (1-dim) 

• Find the conditions for the ground state to be periodic (in any dimension) 

• Prove the devil staircase structure 

• Find the low temperature properties 

• Find the thermodynamic properties of the general class of system consisting of quan- 
tum particles interacting with classical fields. 



A large number of references concerning ltjie Falicov-Kimball model in 1 and 2 dimensions 
can be found in [1] and [2] below and will not be reproduced here, unless quoted in the 
text. 
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